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We perform the careful study of the Gribov copies problem in SU (2) lattice gauge 
theory for maximal direct and maximal indirect center projections. We find that 
this problem is much more severe than it was thought before. The projected string 
tension is not in agreement with the physical string tension. We also show that the 
particle— like objects, nexuses, might be important for the confinement dynamics. 



1 Introduction 

The old idea about the role of the center vortices in confinement phenomena 



has been jjevived recently with t 
invariant □ and gauge dependent 





le use of lattice regularization. Both gauge 
3 approaches are developed. The gauge de- 
pendent studies were done in several center gauges. The center gauge leaves 
intact center group local gauge invariance. It is believed that gauge depen- 
dent P-vortices defined on the lattice plaquettes are able to locate thick gauge 
invariant center vortices and thus provide the essential evidence for the cen- 
ter vortex picture of confinement. So far 3 different center gauges have been 
used in practical computations: the indirect maximal center (IMC) gauge 
the direct maximal center (DMC) gauge □ and the Laplacian center gauge 
The first two suffer from Gribov copies problem. Many results supporting the 
important role of P-vortices are obtained in these two gauges but the problem 
of Gribov copies effects, we are addressing here, has not been studied prop- 
erly. We. also investigate properties of recently introduced new objects called 
nexusesH'Q or center monopolesutl. One can define nexus in SU(N) gauge the- 
ory as a particle-like object formed by N center vortices meeting at the center, 
with the zero (mod N) net flux. We use P-vortices in the center projection to 
define nexuses in SU(2) lattice gauge theory. 

Inn we have found that in DMC gauge the projected string tension is 
much lower than the physical one which contradicts the earlier claims. We also 
confirmed the observation made int3: there are gauge copies which correspond 



"Talk given by V.G.B. at the International Symposium on Quantum Chromodynamics and 
Color Confinement (Confinement 2000), Osaka, March 2000. 
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to higher maxima of the gauge fixing functional F (see below for definition) 
than usually obtained and at the same time these new gauge copies produce 
P-vortices evidently with no center vortex finding ability. These results are 
discussed in section 2. In section 3 we present our results for IMC gauge. 
Results of sections 2 and 3 were obtained with the relaxation - overrelaxation 
(RO) algorithm. In section 4 we present our preliminary results obtained with 
more effective algorithm, simulated annealing (SA). We show that the use of 
this algorit hm permits to obtain higher maxima and thus to solve the puzzle 
imposed in 1121. Moreover SA algorithm gives the_k>west value of the string 
tension. We also discuss the finite volume effects EH for the projected string 
tension. 

2 Direct maximal center gauge 

The DMC gauge is defined i by the maximization of the following functional: 

p V) = (l™^) = w^\ + v > (!) 

with respect to local gauge transformations, U n ,n is the lattice gauge field, 
V is the lattice volume. DMC gauge condition fixes the gauge up to Z{2) 
gauge transformation. The fixed configuration can be decomposed into Z(2) 
and coset parts: U ntll — Z n ^V n . p , where Z n ^ = signTr[/ nAl . Plaquettes 
constructed from Z n>fl field have values ±1. Those of them taking values — 1 
compose the so called P-vortices. P-vortices form closed surfaces in 4D space. 
Some evidence has been collected, that P-vorticesin the center gauges can serve 
to locate gauge invariant center vortices. In ref. □ the projected Wilson loops, 
Wz(2) ? are computed via linking number of the static quarks trajectories and 
P-vortices. It was found that the string tension, <Jz{2) ? obtained from Wz(2) is 
very close to the physical string tension <Jsu(2) ■ This fact has been called center 
dominance. Another important observation was that the density of P-vortices 
scales in agreement with asymptotic scaling qE3. We inspect these statements 
using careful gauge fixing procedure. 

The problem of the DMC gauge fixing is that the functional F(U) has 
many local maxima. This is the analogue of the Gribov problem in continuum 
gauge theoriesllj. We call configurations corresponding to different local max- 
ima Gribov copies. To perform unambiguous computations one must fix the 
gauge completely, i.e. to find the global maximaEJ. It is impossible to do it 
numerically, and we generate a large number N cop of local maxima, calculate 
observables using configuration corresponding to the highest maximum and 
extrapolate results to N cop — > oo limit. The local maxima are produced by ap- 
plying the RO algorithm to random gauge copies of the original configuration. 
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Table 1: The comparison of <7z(2)i a SU(2) an d P for DMC gauge, RO gauge fixing procedure. 



1 y cop 


(T Z(2)/< J SU(2) 


W( fJ S(7(2)a 2 ) 




(3 = 2.3 


(3 = 2.4 


(3 = 2.5 


(3 = 2.3 


(3 = 2.4 


= 2.5 


3 


0.94(2) 


0.93(2) 


0.98(2) 


1.30(1) 


1.51(1) 


1.74(1) 


20 


0.87(2) 


0.80(2) 


0.83(3) 


1.27(1) 


1.42(1) 


1.61(2) 


oo 


0.82(3) 


0.71(3) 


0.71(3) 


1.24(1) 


1.33(2) 


1.49(2) 



Our computations have been performed on lattices L 4 = 12 4 at j3 = 2.3, 2.4 
and L 4 = 16 4 at (3 = 2.5. At (3 = 2.3, 2.4 {(3 = 2.5) we use 100 (50) statistically 
independent gauge field configurations. In Table 1 we show the ratio of string 
tensions <Jz{2)/ a su{2) &z(2) is computed from the Creutz ratio x(J); 3 < 
J < 4 on 12 4 lattice, and 3 < I < 6 on 16 4 lattice. For N cop = 3 (the number 
of gauge copies used ina) <r Z (2) is close to cr SU ( 2 y But it becomes significantly 
lower for large N cop . Thus RO gauge fixing gives the strong dependence of 
<?z(2) on N cop . In the limit N cop — > oo crz(2) is 20-30 % lower than <Jsu(2)- 

In Table 1 we also show the_ratio 2p/asu(2)d 2 (p is the density of P- 
vortices) . As it is claimed in ref. Ej 2p coincides with the dimensionless string 
tension, ugu^a 2 if plaquettes carrying P-vortices are uncorrelated. Our re- 
sults in Table 1 show that the density of P-vortices does not satisfy this relation. 
We have found that for iV cop = 3 n p is in a good agreement with asymptotic 
scaling, as it was observed before But for N cop — > oo p deviates from the 
two loop asymptotic scaling formula and its dependence on [3 becomes similar 
to that of tJ5(7(2)a 2 . 

We also performed computations using the modified gauge fixing procedure 
(LRO) suggested inl!3: every copy has been first fixed to Landau gauge, and 
then the RO algorithm for DMC gauge has been applied. With this procedure 
we found the local maximum higher than local maxima of RO procedure for 
any N cop S [1; 20] and in the limit N cop —* oo (see Table 2). We confirm that 
RO and LRO gauge fixing procedures generate two classes of gauge equivalent 
copies with different properties: <Jz(2) is zero for local maxima produced with 
LRO and nonzero for those produced with RO. The density of P-vortices is 
essentially lower for LRO. Then the idea of the complete gauge fixing based 
on the choice of the global maximum would force us to choose the LRO local 
maxima. Thus the vortex finding property is completely lost. In chapter 4 we 
demonstrate that SA algorithm solves this problem. 



Table 2: < F max > for DMC gauge obtained with various algorithms 





(3 = 2.3 


[3 = 2.4 


(3 = 2.5 


RO{N cop = 


oo) 


0.7552(1) 


0.7764(2) 


0.7943(2) 


LRO(JV cop = 


00 ) 


0.7564(1) 


0.7759(3) 


0.7955(2) 


SA(N cop = 


3) 


0.7588(1) 


0.7770(2) 


0.7970(2) 



2. 1 Nexuses. 

We also investigate the properties of the point like objects, called nexuses. 
On the 4D lattice we have the conserved currents of nexuses, defined af- 
ter the center projection. First we calculate the phase, s/, of the Z(2) link 
variable: Z\ = exp(i7rs;), s; = 0,1. Then we define the plaquette variable 
op = ds mod 2, (erp = 0,1). The nexus current (or center monopole cur- 
rent cl) is then defined as *j — ^S*ap. These currents live on the surface of 
the P-vortex (on the dual 4D lattice) and P-vortex flux goes through positive 
and negative nexuses in alternate order. The important characteristic of the 
cluster of currents is the condensate, C, definedli3 as the percolation probabil- 
ity. As it is shown in ref. a the condensate C of the nexus currents is the order 
parameter for the confinement - deconfinement phase transition. We found 
that C is nonzero for the gauge copies obtained via RO procedure (when the 
projected Wilson loops have the area law). On the other hand C is zero (in the 
thermodynamic limit L — > oo) for gauge copies obtained using LRO procedure 
(when the projected Wilson looos have no area law). It is interesting that for 
RO procedure C seems to scale a as the physical quantity with the dimension 
(mass)' 1 . This is illustrated in Fig.l, where we plot the /3-dependence of the 
ratio C / (a su (2)a 2 ) 2 ■ Thus these new objects might be important degrees of 
freedom for the description of the nonperturbative effects. 

3 Indirect maximal center gauge 

The gauge fixing condition for IMC gauge consists of two steps : 
Maximal abelian gauge (SU (2) — > U(l)) is fixed by solving maximization prob- 
lem : 

max{iq(C/»)}, g G SU(2)/U(1); F 1 (U) = -J-^Tr (U ntlt a 3 ^a a ) 
and U(l) field is extracted: U n ^ = V n ^u n ^ , u n _ p £ U(\). 
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Figure 1: The /3-dependence of the ratio of the nexus condensate, C, to the SU(2) string 
tension in lattice units. 



Then Z(2) gauge is fixed by maximizing 

max{F 2 K)}; w 6 U(l)/Z(2), F 2 (u) = ^ £ (Rcu n ^) 2 . 

The Z(2) field is defined as: u n _^ ~ Z n ^v n ^. 

In this section we discuss our results obtained with two gauge fixing pro- 
cedures: 

- RO procedure: RO algorithm at both steps; 

- LRO procedure: LRO algorithm at the 1st step and RO algorithm at the 
2nd one; 

Our results for the ratios o'z(2)/' J su{2) an d 2p/csu{2) a2 obtained with RO 
procedure are presented in the Table 3. These results arc in good qualitative 
and even quantitative agreement with the results obtained in DMC gauge. The 
LRO procedure (as in DMC gauge case) corresponds to uz(2) = 0. Another 
important observation is that at the first stage (MA gauge) LRO procedure 
local maxima average < Fi^ max > is higher than those of RO procedure (see 
Table 4). Thus the problem discussed in the previous section is also relevant 
for IMC gauge. 
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Table 3: The comparison of <tz(2)i a SU(2) an d P f° r IMC gauge, RO gauge fixing procedure. 





fJ Z(2)/c r SC/(2) 




■*Y cop 


/3 = 2.3 


(3 = 2.4 


= 2.5 


(3 = 2.3 


(3 = 2.4 


= 2.5 


(1,1) 


1.03(1) 


1.06(1) 


1.09(1) 


1.51(1) 


1.84(1) 


2.17(1) 


(20,10) 


0.89(3) 


0.81(2) 


0.86(2) 


1.44(1) 


1.69(1) 


1.97(2) 


(oo,oo) 


0.81(4) 


0.69(3) 


0.72(3) 


1.40(1) 


1.63(1) 


1.83(2) 



Table 4: < Fi yJnax > for IMC gauge obtained with various algorithms 





f3 = 2.3 


f3 = 2.4 


[3 = 2.5 


RO {N cop = oo) 


0.7132(1) 


0.7324(3) 


0.7509(2) 


LRO(A cop = oo) 


0.7134(1) 


0.7331(2) 




SA(A cop = oo) 




0.7337(2) 





4 Results of the SA algorithm . 

The simulated annealing algorithm has been applied to the maximal abelian 
gauge fixing intZUla and its advantages in reducing the bias due to Gribov copies 
effects has been demonstrated. In this section we describe our preliminary 
results in fixing for both DMC and IMC gauges obtained with SA algorithm. 
For IMC gauge SA algorithm has been applied at the first step only, while RO 
algorithm has been applied at the 2nd step. The details of the SA algorithm 
for DMC gauge will be explained elsewhere. 

Our first observation is that with SA algorithm it is possible to reach higher 
local maxima < F max > (for DMC gauge) and < Fi, max > (at the first step 
of IMC gauge - MA gauge) than those reached with the LRO procedure. This 
can be seen from the Tables 2 and 4. Our second observation is that results 
obtained with these local maxima are qualitatively similar to those obtained 
with RO procedure, i.e. uz(2) computed on these configurations is nonzero. 
The values of <Jz(2) and p are lower than those obtained with RO algorithm. 

For SA algorithm applied to DMC gauge we generated only 3 local maxima 
per configuration, since this algorithm is rather time consuming. Thus we 
are not able to make extrapolation N cop — > oo. Nevertheless the example 
N cop = 3 is instructive. Our results are = 0.69(2), 0.72(2), 0.49(2) at 

(3 = 2.3, 2.4, 2.5 correspondingly. Comparing with results in Table 1 one can 
see that already for N cop — 3 SA values for az2 are lower than values obtained 
with RO procedure in N cop — oo limit. For IMC gauge with SA algorithm 
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at the first step we get = 0.73(4) at /3 = 2.4 (computations have been 

made at this (3 only) . This is in agreement with the result in Table 3 obtained 
with RO procedure. It is still possible that employing SA algorithm also at the 
second step of IMC gauge fixing will bring further decreasing of the projected 
string tension. 

IntLil it has been argued that low values for <tz2 reported inlj can be due 
to finite volume effects which spoil vortex finding property even on the lattices 
where finite volume effects are not visible in the gauge invariant observables. 
We made some computations for DMC gauge with RO procedure for lattices 
10 4 to 20 4 and found out results in qualitative agreement with fig. 3 of ref. E3J. 
But the problem is still not settled since increasing of az2 with lattice volume 
can be due to lower value of local maxima found on larger lattices where finding 
higher maxima becomes too costly^ while there is clear anticorrelation between 
az2 and the value of < F max > tl. Moreover our results with SA algorithm 
indicate that these finite volume effects must be much more pronounced if one 
takes the highest maxima available. We made computation in DMC gauge on 
L = 16 lattice at (3 — 2.4. We apply the SA algorithm with more updating 
sweeps than it was done for L = 12 computations at the same j3. We obtained 
higher value < F max >= 0.7784(1) and lower value = °- 64 ( 2 ) than 

corresponding values given in Table 1. The question of the finite volume effects 
for IMC gauge has not been studied so far. 

Conclusions 

We conclude that DMC and IMC gauges suffer from strong Gribov copies 
effects and careful gauge fixing is necessary to make the bias caused by these 
effects reasonably small. This procedure is rather costly. Another problem of 
DMC and IMC gauges is that (Jz(2) is too low. It is not clear whether relation 
f7 z(2) = < J su(2) can be valid at large enough lattices. Even if this is the case, 
the size of these lattices is enormous. The alternative gauge, the Laplacian 
center gauged is then more favorable. 
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